Location: CFTR @ c7f30a71ce7f / CFTR - BG.cellml

Author:
Glucose Uptake <nafs080@aucklanduni.ac.nz>
Date:
2017-09-08 10:09:49+12:00
Desc:
Initial model of CFTR
Permanent Source URI:
https://staging.physiomeproject.org/workspace/4c0/rawfile/c7f30a71ce7f2662d885269bb316241003f87b65/CFTR - BG.cellml

<?xml version='1.0'?>
<model name="CFTR_BG" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#" xmlns:xlink="http://www.w3.org/1999/xlink">
    <import xlink:href="Units.cellml">
        <units name="per_mol" units_ref="per_mol"/>
        <units name="J_per_mol" units_ref="J_per_mol"/>
        <units name="mol_per_s" units_ref="mol_per_s"/>
        <units name="Js_per_C2" units_ref="Js_per_C2"/>
        <units name="C_per_mol" units_ref="C_per_mol"/>
    </import>
    <import xlink:href="Parameters.cellml">
        <component component_ref="param" name="param"/>
    </import>
    <component name="CFTR">
        <variable name="t" public_interface="out" units="second"/>
        <!-- Parameters-->
        <variable name="RT" public_interface="in" units="J_per_mol"/>
        <variable name="R" public_interface="in" units="Js_per_C2"/>
        <variable name="U_E" public_interface="in" units="J_per_mol"/>
        <!--        var K_Gl_o: per_mol {pub: in};
        var K_Gl_i: per_mol {pub: in};-->
        <!--var K_HCO3_i: per_mol {pub: in};
var K_HCO3: per_mol {pub: in};-->
        <!--        var K_E_o: per_mol {pub: in};
        var K_EGl_o: per_mol {pub: in};
        var K_EGl_i: per_mol {pub: in};
        var K_E_i: per_mol {pub: in};-->
        <!--var K_EHCO3_i: per_mol {pub: in};
var K_EHCO3_o: per_mol {pub: in};-->
        <!--        var K_Re1: mol_per_s {pub: in};
        var K_Re2: mol_per_s {pub: in};
        var K_Re3: mol_per_s {pub: in};
        var K_Re4: mol_per_s {pub: in};-->
        <!--var K_Re5: mol_per_s {pub: in};
var K_Re6: mol_per_s {pub: in};-->
        <!-- State variables-->
        <variable initial_value="1.0" name="q_Cl_i" units="mole"/>
        <variable initial_value="0.1" name="q_Cl_o" units="mole"/>
        <!--        var q_HCO3_i: mole {init: 0.0};
        var q_HCO3_o: mole {init: 1.0};-->
        <!--        var q_E_i: mole {init: 1.0};
        var q_EGl_i: mole {init: 0.0};
        var q_EGl_o: mole {init: 0.0};
        var q_E_o: mole {init: 0.0};-->
        <!--        var q_EHCO3_i: mole {init: 0.0};
        var q_EHCO3_o: mole {init: 0.0};-->
        <variable name="v_Re1" public_interface="out" units="mol_per_s"/>
        <!--        var v_Re2: mol_per_s {pub: out};
        var v_Re3: mol_per_s {pub: out};
        var v_Re4: mol_per_s {pub: out};-->
        <!--        var v_Re5: mol_per_s {pub: out};
        var v_Re6: mol_per_s {pub: out};-->
        <!-- Conservation laws-->
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>q_Cl_i</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v_Re1</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>q_Cl_o</ci>
                </apply>
                <ci>v_Re1</ci>
            </apply>
            <!--        ode(q_HCO3_i, t) = -v_Re4;
        ode(q_HCO3_o, t) = v_Re6;-->
            <!--        ode(q_E_i, t) = v_Re3-v_Re1;
        ode(q_EGl_i, t) = v_Re1-v_Re4;
        ode(q_EGl_o, t) = v_Re4-v_Re2;
        ode(q_E_o, t) = v_Re2-v_Re3;-->
            <!--        ode(q_EHCO3_i, t) = v_Re4-v_Re5;
        ode(q_EHCO3_o, t) = v_Re5-v_Re6;-->
            <!-- Constitutive relations-->
            <apply>
                <eq/>
                <ci>v_Re1</ci>
                <apply>
                    <times/>
                    <ci>R</ci>
                    <apply>
                        <minus/>
                        <ci>U_E</ci>
                        <apply>
                            <times/>
                            <ci>RT</ci>
                            <apply>
                                <ln/>
                                <apply>
                                    <divide/>
                                    <ci>q_Cl_i</ci>
                                    <ci>q_Cl_o</ci>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
        </math>
        <!--        v_Re2 = K_Re2*(K_EGl_o*q_EGl_o-K_Gl_o*q_Gl_o*K_E_o*q_E_o);
        v_Re3 = K_Re3*(K_E_o*q_E_o-K_E_i*q_E_i);
        v_Re4 = K_Re4*(K_EGl_i*q_EGl_i-K_EGl_o*q_EGl_o);-->
        <!--        v_Re5 = K_Re5*(K_EHCO3_i*q_EHCO3_i-K_EHCO3_o*q_EHCO3_o);
        v_Re6 = K_Re6*(K_EHCO3_o*q_EHCO3_o-K_HCO3_o*q_HCO3_o*K_E_o*q_E_o);-->
    </component>
    <!-- Mapping to the Parameters module-->
    <connection>
        <map_components component_1="param" component_2="CFTR"/>
        <map_variables variable_1="RT" variable_2="RT"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="U_E" variable_2="U_E"/>
        <!--        vars K_Gl_o and K_Gl_o;
        vars K_Gl_i and K_Gl_i;-->
        <!--        vars K_HCO3_i and K_HCO3_i;
        vars K_HCO3_o and K_HCO3_o;-->
        <!--        vars K_E_o and K_E_o;
        vars K_EGl_o and K_EGl_o;
        vars K_EGl_i and K_EGl_i;
        vars K_E_i and K_E_i;-->
        <!--        vars K_EHCO3_i and K_EHCO3_i;
        vars K_EHCO3_o and K_EHCO3_o;-->
        <!--        vars K_Re1 and K_Re1;
        vars K_Re2 and K_Re2;
        vars K_Re3 and K_Re3;
        vars K_Re4 and K_Re4;-->
        <!--        vars K_Re5 and K_Re5;
        vars K_Re6 and K_Re6;-->
    </connection>
</model>